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The Vibrionaceae (Vibrio) are a ubiquitous group of metabolically flexible marine bacteria 
that play important roles in biogeochemical cycling in the ocean. Despite this versatility, 
little is known about Vibrio diversity and abundances in upwelling regions. The seasonal 
dynamics of Vibrio populations was examined by analysis of 16S rRNA genes in Monterey 
Bay (MB), California from April 2006— April 2008 at two long term monitoring stations, 
C1 and M2. Vibrio phylotypes within MB were diverse, with subpopulations clustering 
with several different cultured representatives including Allivibrio spp., Vibrio penaecida, 
and Vibrio splendidus as well as with many unidentified marine environmental bacterial 
16S rRNA gene sequences. Total Vibrio population abundances, as well as abundances 
of a Vibrio sp. subpopulation (MBAY Vib7) and an Allivibrio sp. subpopulation (MBAY Vib4) 
were examined in the context of environmental parameters from mooring station and CTD 
cast data. Total Vibrio populations showed some seasonal variability but greater variability 
was observed within the two subpopulations. MBAY Vib4 was negatively associated with 
MB upwelling indices and positively correlated with oceanic season conditions, when 
upwelling winds relax and warmer surface waters are present in MB. MBAY Vib7 was 
also negatively associated with upwelling indices and represented a deeper Vibrio sp. 
population. Correlation patterns suggest that larger oceanographic conditions affect the 
dynamics of the populations in MB, rather than specific environmental factors. This study 
is the first to target and describe the diversity and dynamics of these natural populations 
in MB and demonstrates that these populations shift seasonally within the region. 
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INTRODUCTION 

The Vibrionaceae (Vibrio) are a group of physiologically- flexible 
marine bacteria that are ubiquitous in ocean waters and have 
been identified in most marine ecosystems (Wietz et al., 2010). 
Vibrio have the distinctive ability to break down and utilize many 
carbon, nitrogen, and phosphorus substrates (McDougald and 
Kjelleberg, 2006; Thompson and Polz, 2006; Dryselius et al., 
2008; Lai et al, 2009; Salter et al, 2009) and their production of 
the external enzymes chitinase and laminarinase provide access 
to abundant nutrients that are unavailable to other organisms 
(Svitil et al, 1997; Riemann et al., 2000; Ansede et al, 2001; 
Alderkamp et al., 2007; Murray et al, 2007). In addition to 
their diverse metabolic capabilities, Vibrio species have devel- 
oped adaptive responses to starvation and environmental stress 
which include conversion to an ultramicrobacterial morphology 
(<0.4(jim diameter) (Denner et al., 2002) and retention of a 
high concentration of rRNA as a "stimulation ready" response 
that may contribute to this group's rapid growth potential after 
nutrient influxes (Eilers et al., 2000). Despite these capabili- 
ties, there have been few studies in oceanic upwelling regions, 
where metabolic flexibility, rapid nutrient response, and highly 
developed stress protection mechanisms should prove to be eco- 
logically advantageous during changing environmental regimes. 

Monterey Bay (MB) is an open embayment along the 
California coast that is distinguished by a near-shore deep-water 



canyon where periodic upwelling events sustain diverse sea life. 
Circulation within MB is variable and influenced by recently 
upwelled waters as well as offshore waters from the California 
and Davidson Currents that enter MB during relaxation events 
(Rosenfeld et al., 1994). Three hydrographic seasons have been 
defined within MB — an upwelling or "cold water phase" that 
usually occurs from mid-February through August, an oceanic 
period or "warm water phase" from mid-August through mid- 
October and a Davidson Current or "low thermal gradient" 
phase between mid-November and mid-February (Skogsberg and 
Phelps, 1946; Breaker and Broenkow, 1994). The periodic influx 
of nutrient-enriched waters into MB via upwelling results in a 
highly productive ecosystem presenting a unique environment for 
the study of Vibrio population dynamics. 

Vibrio species are common isolates from MB waters but most 
studies have concentrated on pathogenic representatives (Kenyon 
et al, 1984; Kaysner et al, 1987; Miller et al, 2006). One phyloge- 
netic screening study found that Vz'fon'o-specific sequences made 
up 3.1% of the genetic information within a FOSMID library 
developed from a 100 m deep sample (Suzuki et al, 2004). Vibrio 
populations are estimated to average only about 1% of marine 
bacterial populations worldwide (Thompson and Polz, 2006) but 
it has also been suggested that differences in vertical distributions 
may be more significant than differences in geographic distri- 
bution for Vibrio populations, and thus, may make up a more 
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significant portion of subsurface bacterial populations (Simidu 
and Tsukamoto, 1985). 

This study was designed to determine the seasonal variability 
of Vibrio populations in MB by examining total Vibrio pop- 
ulation dynamics as well as two subpopulations, a Vibrio sp. 
subpopulation (MBAY Vib7) and an Allivibrio sp. subpopulation 
(MBAY Vib4), and to analyze how population diversity changes 
with upwelling, season, and environmental factors. MB has a sys- 
tem of long term monitoring stations (Pennington and Chavez, 
2000) and two of the stations within this network, CI and M2 
(Figure 1), were chosen as sample sites for this study. CI is a 
coastal site that is influenced by aged upwelled waters and sum- 
mer relaxation events. M2 is an outer Bay site that is influenced by 
upwelled waters that mix with California Current and Davidson 
Current waters. The examination of Vibrio population dynamics 
at these two different locations can provide insight into the effect 
of upwelling characteristics on Vibrio abundances. 

MATERIALS AND METHODS 
SAMPLE COLLECTION 

Seawater samples were collected from MB at stations CI 
(36.797N; 121.847 W; Figure 1) and M2 (36.697N; 122.378 W; 
Figure 1) during cruises of the MB Time Series project under- 
taken by the Biological Oceanography Group (BOG) at the 
Monterey Bay Aquarium Research Institute (MBARI). Samples 
were collected from April 2006 through April 2008 on a peri- 
odic basis spanning 19 sampling dates. Samples were collected 
using a SeaBird 911 CTD rosette equipped with physical sen- 
sors described by Pennington and Chavez (2000). A total of 82 
samples were collected at station CI at 5 m (15 samples), 10 m 
(17 samples), 20 m (18 samples), 30 m (17 samples), and 200 m 
(15 samples). At station M2 a total of 86 samples were col- 
lected at 5 m (12 samples), 10 m (14 samples), 20 m (14 samples), 
40 m (16 samples), 100 m (16 samples), and 200 m (16 samples). 
For each sample 1-2 liters of seawater were filtered using gentle 




122.4 W 122.2 W 122 W 121.8 W 



FIGURE 1 | Map of sample sites C1 and M2 in Monterey Bay, California. 

Each individual blue dot represents a sampling event and the black circle is 
the identified coordinate of the mooring station at M2. 



peristaltic pumping through sequential in-line 25 mm diame- 
ter 10 |xm pore-size PE filters (GE Osmonics, Minnetonka, MN, 
USA) and 25 mm diameter 0.2 |i,m pore-size Supor 200 mem- 
brane filters (Pall Corporation, Port Washington, NY, USA). The 
filters were transferred to 1.5 mL polypropylene microcentrifuge 
tubes containing 0.2 g of 0. 1 mm and 0.5 mm diameter autoclaved 
glass beads (BioSpec Products, Bartlesville, OK, USA). Samples 
were flash frozen in liquid nitrogen before transfer to a — 80°C 
freezer onshore for storage until nucleic acids were extracted. 

ENVIRONMENTAL DATA 

Upwelling indices (UI) of both monthly and daily averaged 
upwelling conditions were obtained for 36°N 122°W from the 
Pacific Fisheries Environmental Laboratory (PFEL, Pacific Grove) 
for determination of upwelling patterns (http://www.pfeg.noaa. 
gov/products/PFEL/modeled/indices/PFELindices.html) . Units 
are given as m 3 s _1 100 m of coastline- 1 as the average amount 
of water upwelled through the bottom of the Ekman layer each 
second along each 100 m of coastline on a scale of about 200 
miles. Surface data for nitrate, chlorophyll a, and temperature 
was obtained from the M2 mooring from MBARI LOBOviz 
(http://www.mbari.org/lobo/loboviz.htm). The BOG at MBARI 
provided environmental data from CTD measurements as well as 
surface ( < 10 m) phytoplankton concentrations, which were ana- 
lyzed as described by Chavez et al. (1990). Specific phytoplankton 
groups (Synechococcus, total diatom, total dinoflagellate, and total 
phytoplankton) were chosen for analysis to assess the influence 
of different phytoplankton regimes on the variability of Vibrio 
groups in the water column. 

DNA EXTRACTION 

DNA was extracted from the 0.2 \im pore-size filters using the 
Qiagen DNeasy Plant Kit (Hilden, Germany) as described by 
Moisander et al. (2008) with modifications to improve DNA 
recovery. Cells were lysed by a triplicate run of freeze-thaw steps in 
liquid nitrogen followed by a 65° C water bath. Cells were further 
disrupted by a 2-min agitation of bead beating (Mini bead beater 
96, BioSpec Products, Bartlesville, OK, USA) and DNA yield was 
increased with the addition of 0.45 |xL of Proteinase K and incu- 
bation at 55° C for 1 h. AE buffer was used as the elution medium. 
Two duplicate elutions of 25 \iL each were combined for a final 
elution volume of 50 u,L. Extracts were stored at — 20° C until use. 

AMPLIFICATION, CLONING AND SEQUENCING OF l/ifirio-SPECIFIC 16S 
rRNA GENE 

Samples were selected from 2007 to 2008 from 5 to 200 m 
depths at both station CI and M2 to construct clone libraries 
of amplified 16S rRNA gene sequences from the MB. A total 
of 65 samples were processed for inclusion in the libraries. 
The PCR amplification utilized a universal forward primer (27F, 
Table 1) and a Vibrio specific reverse primer (680R, Table 1) 
and followed the two-phase PCR amplification technique out- 
lined by Thompson et al. (2004b). Samples were amplified using 
Invitrogen Taq polymerase (Carlsbad, California), on a BioRad 
thermal cycler (Hercules, California). The PCR products were gel 
purified and cloned into pGEM-T vectors (Promega, Madison 
WI) using manufacturer's guidelines. Ten to fifteen clones were 
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Table 1 | Primers and probes for PCR and qPCR (written 5'->3'). 



16S rRNA GENE (£. coli POSITION 27-680) 


27F 


AGAGTTTG ATC MTGG CTCAG 


680 R 


G AAATTCTAC C C C CCTCTAC AG 




Forward (587F) 


GGCGTAAAGCGCATGCAGGT 


Reverse (680R) 


G AAATTCTAC C C C CCTCTAC AG 




Forward (Vib4F) 


GGGAGGAAGGTGTTGTCGTTAA 


Reverse (Vib4R) 


ACGGAGTTAGCCGGTGCTT 


Probe (Vib4P) 


[FAMIAGCGGCAGCATTTGACGTTACCCACITAMRA] 




Forward (Vib7F) 
Reverse (Vib7R) 
Probe (Vib7P) 


G CAGTGAG GAAG G GTGTGTAGT 

GAG TTAG CC G GTGCTTCTTCTG 

[FAM]ATAGCTGCATACTTTGACGTTAGCTG[TAMRAl 



The 16S rRNA gene and Total Vibrio primers are from Thompson et al. (2004a, b). 
The MBAY Vib4 and MB AY Vib7 primer and FAM/TAMRA labeled probe sets 
were developed based on sequences from the Monterey Bay 16S rRNA gene 
clone libraries. 



chosen from each sample, for a total of 911 clones, and pre- 
pared for sequencing with the Montage Plasmid Miniprep kit 
(Millipore, Billerica, MA) according to manufacturer's proto- 
cols. Cloned inserts were sequenced at the UC Berkeley DNA 
Sequencing Facility using T7 primers and analyzed on an Applied 
Biosystems 3730x1 DNA Analyzer. Sequences were aligned and 
compared to published sequences using the Ribosomal Database 
Project (RDP) on-line interface (Cole et al, 2007, 2009) and 
were quality checked for chimeras using the RDP Chimera check 
program and Bellerophon (Huber et al., 2004). Phylogenetic anal- 
ysis of 827 sequences was conducted in ARB (Ludwig et al., 
2004) and neighbor-joining phylogenetic trees were constructed 
using the Jukes-Cantor correction. The distance matrix derived 
from the neighbor joining analysis was used in DOTUR for 
assignment of operational taxonomic units (OTUs) (Schloss and 
Handelsman, 2005). All sequences were submitted to the National 
Center for Biotechnology Information (NCBI) GenBank database 
as accession numbers KF941543-KF942369. 

QUANTITATIVE PCR (qPCR) 

Total Vibrio abundances were quantified using a SYBR Green 
qPCR method utilizing the Vibrio specific primers 567F and 
680R (Table 1) from Thompson et al. (2004b) to ensure a broad 
spectrum of Vibrio species were included in the analysis. DNA 
template (2 |xl) was added to 12.5 u,l Applied Biosystems SYBR® 
Green PCR Master Mix, 8.25 |xl of 5 kD filtered water, 0.25 (xl of 
bovine serum albumin (BSA) and 1 u,l of each primer for a total 
volume of 25 |xl. The cycling conditions were 50° C for 2min, 
95°C for 10 min followed by 40 cycles of 95°C for 15 s and 64°C 
for 1 min. Each run was followed by a dissociation step (95°C 
for 15 s and 60°C for 1 min 95°C for 15 s and 60°C for 15 s) to 
determine a melt curve for analysis of specificity. 

Two qPCR assays utilizing TaqMan® chemistry were designed 
to target two Vibrio subpopulations within MB based on 
sequences of interest identified from the MB 16S rRNA clone 
library. The MBAY Vib4 group was closely associated with 



Alllivibrio sp. and the MBAY Vib7 group was associated with 
Vibrio penaecida. Primers and probes were designed with Primer 
Express (Table 1) (ABI, Foster City, CA) and synthesized by 
Sigma-Genosys (Woodlands, TX). The reactions for the sub- 
population assays contained 2 \i\ of sample DNA, 12.5 u,l Applied 
Biosystems TaqMan® Universal PCR Master Mix, 8 u.1 of 5 kD 
filtered water, 1 \i\ each of the forward and reverse primer and 
0.5 ul of the probe. Samples were assessed using a TaqMan® 
qPCR method with cycling conditions of 50° C for 2 min, 95° C 
for 10 min followed by 45 cycles of 95°C for 15 s and 60°C for 
1 min. The qPCR primer sets were used to analyze the full range 
of samples at both stations. 

Both qPCR assays compared Cj values to standard curves 
(equivalent to 10 1 — 10 s gene copies per reaction) derived from 
group specific environmental isolate clones from the MB clone 
library linearized with Nadl. An internal control at a concentra- 
tion of 10 4 gene copies of linearized standard assessed the runs for 
inhibition. If detected, DNA samples were diluted 1:10 and ampli- 
fication was repeated. All qPCR runs were performed on a 7500 
Applied Biosystems Real-Time PCR instrument. Final concentra- 
tions are reported as gene copies l - 1 of seawater sampled. 

STATISTICAL ANALYSIS 

Statistical analysis was performed in JMP version 9 (SAS Institute, 
North Carolina). Data was natural log transformed [In (x+ 1)] 
to scale variables for graphical interpretation and adjust for 
normalcy. A Principal Components Analysis (PCA) bi-plot was 
utilized to visualize relationships between Vibrio populations and 
physical and biological variables and Spearman's rank analysis 
(p s ) was used to define statistically significant relationships. 

RESULTS 

ENVIRONMENTAL CONDITIONS IN MB 

Water temperatures from the sample sites ranged from 7.7 to 
15.4°C (Table 2). Salinity ranged from 32.7 to 34.4 ppt and 
chlorophyll a values ranged from 0 to 18.45 (xg liter -1 (Table 2). 
The Monthly Upwelling Index (MUI) ranged from 46 to 237 m 3 
s~ 1 1 00 m of coastline - 1 . The daily Upwelling Index was analyzed 
against surface conditions at M2 (Table 3A). Results were consis- 
tent with expected trends for upwelling events with a negative cor- 
relation to temperature (p s = —0.47, p = 0.002) and a positive 
correlation to both nitrate (p s = 0.22, p = 0.035) and chlorophyll 
a (p s = 0.28, p = 0.009). Principal Component Analysis asso- 
ciated total diatoms, total dinoflagellates and all phytoplankton 
with the upwelling season and the monthly UI. Synechococcus was 
more closely associated with station M2 (Figure 2). 

16S rRNA GENE PHYLOGENETIC ANALYSIS 

Multiple distinct 16S rRNA-based phylotypes were recovered 
from the MB clone libraries (Figure 3). The primer sets uti- 
lized to develop the MB clone libraries targeted Vibrio sp. 
strains but 67 sequences within the library clustered with the 
closely related Photobacterium sp. There were also a significant 
number of sequences that clustered with unidentified Vibrio 
strains and uncultured marine bacteria. Vibrio splendidus was 
defined in several distinct clades within the phlogenetic tree. 
Other MB sequences clustered with V. pomeroyi and V. lentus; 
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Table 2 | CTD measurement and Vibrio concentration ranges over the study period for all stations and depths. 



Station 


Depth 


Temperature 


Salinity 


chlorophyll a 


Total Vibrio 


MBAY Vib4 


MBAY Vib7 


<°C) 


(PPt) 




(gene copies I 1 seawater) 




A II 
All 


All 


"7 "7 1C y| 


QO "7 OA A 

Jz. /— J4.4 


U.UU— I o.4o 


1.80 x 10 3 -3.72 x 10 6 


nd-1.95 x 10 6 


na-o. I o x iu 


C1 


5 


10.7-14.7 


33.4-33.9 


0.76-18.45 


2.14 x 10 4 -3.72 x 10 6 


nd-1.95 x 10 6 


nd-1.07 x 10 4 




10 


10.5-14.3 


32.7-34.0 


0.88-6.61 


2.20 x 10 4 -2.28 x 10 6 


nd-4.03 x 10 5 


nd-2.22 x 10 4 




20 


9.4-13.7 


32.9-33.9 


0.01-9.17 


2.05 x 10 3 -2.83 x 10 6 


nd-1.54 x 10 6 


nd-1.69 x 10 4 




30 


8.9-13.4 


33.2-34.0 


0.00-2.87 


5.08 x 10 3 -2.57 x 10 6 


nd-1.02 x 10 5 


nd-5.05 x 10 4 




200 


7.7-10.0 


33.9-34.2 


0.00-0.34 


1.05 x 10 5 -1 .90 x 10 6 


nd-4.83 x 10 5 


nd-3.92 x 10 4 



M2 5 10.7-15.4 33.3-34.0 0.39-8.59 3.42 x 10 3 -3.13 x 10 6 nd-3.99 x 10 4 nd-9.35 x 10 4 

10 10.5-15.1 33.3-33.9 0.60-9.84 2.23 x 10 3 -1 .01 x 10 6 nd-2.95 x 10 4 nd-5.50 x 10 4 

20 10.5-14.6 33.2-33.7 0.40-5.18 1 .80 x 10 3 -2.58 x 10 5 nd-2.68 x 10 4 nd-6.58 x 10 4 

40 9.7-13.3 33.0-33.8 0.05-1.38 2.18 x 10 4 -4.77 x 10 5 nd-5.78 x 10 4 nd-1.09x10 5 

100 8.6-10.6 33.7-34.0 0.00-0.34 7.18 x 10 3 -7.13 x 10 s nd-3.02 x 10 4 nd-5.18x10 5 

200 7.9-9.4 34.0-34.4 0.00-0.17 2.29 x 10 4 -6.92 x 10 5 nd-1.63 x 10 4 nd-1.07 x 10 s 



nd, none detected. 



V. pectenicida and V. rumoiensis; and V. aesturianus. One clade 
included sequences that clustered with V. agarivorans, V. hispani- 
cus, V. haliticoli, V. ezurae, V. fortis, V. proteolytics, V. sinaloensis, 
V. nigripulchritudo, V. harveyi, V. natriegens, and V. alginolyticus 
(Assorted Vibrio in Figure 3). Quantitative PCR (qPCR) assays 
were developed for two specific phylotypes recovered from the 
MB clone libraries. The MBAY Vib4 group (Figure 3B) clustered 
with Allivibrio sp., which includes A. salmonicida, A. fischeri, A. 
wodanis, and A. logei. MBAY Vib4 included sequences obtained 
from samples collected from both stations CI and M2. Within the 
clone library the MBAY Vib7 group included sequences that were 
only obtained at station M2 and clustered with Vibrio penaecida 
(Figure 3C). 

Rarefaction curves were determined based on the neighbor- 
joining distance matrix of sequences within the MB clone 
libraries. Of the 827 sequences included in the total sequences 
analysis, 158 OTUs were defined at 97% identity, 103 OTUs at 
95% identity, 61 OTUs at 93% identity, and 30 OTUs at 90% 
identity. Only at 93 and 90% identity did the rarefaction curves 
begin to reach an asymptote. Analysis revealed that there was no 
significant difference in species richness between station CI and 
M2 but that both sampling efforts did not attain full coverage to 
estimate total population OTUs. Rarefaction analysis of the 327 
sequences from station CI defined 97 OTUs at 97% identity, 69 
OTUs at 95% identity, 45 OTUs at 93% identity, and 23 OTUs 
at 90% identity. Rarefaction analysis of the 500 sequences from 
station M2 defined 98 OTUs at 97% identity, 66 OTUs at 95% 
identity, 41 OTUs at 93% identity, and 21 OTUs at 90% identity. 

QUANTIFICATION OF Vibrio POPULATION ABUNDANCES 

There were observed differences in the Vibrio populations at both 
of the stations, with the MBAY Vib4 group having higher peak 
abundances at station CI (Figure 4). The MBAY Vib4 group 
ranged in concentration from undetectable levels at all depths 
and stations up to 1.95 x 10 6 gene copies liter 1 seawater at 5 m 
at station CI (Table 2). MBAY Vib4 was positively correlated to 
temperature (p s = 0.32, p = 0.002), oceanic season (p s = 0.41, 



p = 0.002), Synechococcus (p s = 0.20, p = 0.02), and total Vibrio 
abundance (p s = 0.55, p < 0.001) (Table 3). This group was 
negatively correlated to salinity (p s = —0.22, p = 0.02), chloro- 
phyll a (p s = -0.16, p = 0.035), monthly UI (p s = -0.61, p < 
0.001), surface phytoplankton ( p s = -0.44,p < 0.001), and total 
diatoms (p s = -0.55, p < 0.001) (Table 3). 

The MBAY Vib7 population was not as abundant as MBAY 
Vib4 (Figure 4), with concentrations ranging from undetectable 
to 5.18 x 10 5 gene copies liter -1 seawater (Table 2). MBAY 
Vib7 was positively correlated with depth (p s = 41, p < 0.001) 
(Table 3), with the highest detected concentration at 100 m 
at station M2 (Table 2). MBAY Vib7 populations were posi- 
tively associated with station M2 (p s = 0.53, p < 0.001), sur- 
face Synechococcus populations (p s = 0.41, p < 0.001), and the 
Davidson season (p s = 0.15, p = 0.0486) (Table 3). Temperature 
(p s = -0.30, p = 0.0005), chlorophyll a (p s = -0.47, p < 0.001) 
and monthly UI (p s = — 0.23,p = 0.012) were all negatively asso- 
ciated with MBAY Vib7 (Table 3). MBAY Vib7 was identified in 
samples from both station CI and M2 but at lower concentrations 
and at deeper depths at station CI. 

Total Vibrio populations showed some variability in abun- 
dance but were consistently present throughout the study at all 
depths and stations (Figure 4). Total Vibrio population abun- 
dances ranged from 1.80 x 10 3 to 3.72 x 10 6 gene copies liter -1 
seawater (Table 2). The abundance of total Vibrio populations 
was positively correlated to MBAY Vib4 (p s = 0.55, p < 0.001), 
and the oceanic season ( p s = 0.3 l,p < 0.001) and negatively cor- 
related to station number (p s = —0.32, p < 0.001) the monthly 
UI (p s = -0.26, p = 0.0039) and upwelling season (p s = -0.21, 
p = 0.0067) (Table 3). No other parameters examined in this 
study showed significant correlation with the abundance of the 
total Vibrio population. 

DISCUSSION 

Genetic diversity within Vibrio populations can be up to 7% 
when examining the entire 16S rRNA gene (Dorsch et al., 1992; 
Kitatsukamoto et al., 1993). This variability is mostly located 
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FIGURE 2 | Principal Component Analysis of physical and biological 
parameters. MBAY Vib4 (O) MBAY Vib7 (A) from stations C1 (red) and M2 
(blue). 



Table 3 | Spearman's rho (p s ) statistical analysis of physical and 
biological parameters. 



(A) UPWELLING PARAMETERS 








Temperature 


Chlorophyll a 


Nitrate 


Ul 


-0.47 


0.28 


0.22 


Nitrate 


-0.53 


0.15 




Chlorophyll a 


-0.20 










lotal vibrio 


IVIdAY VID4 


HAD AV \/:U"7 

MBAY VlD/ 


Station number 


-0.32 


-0.16 


0.53 


Depth 


0.15 


0.00 


0.41 


Temperature 


-0.03 


0.32 


-0.30 


Salinity 


0.05 


-0.22 


0.12 


Chlorophyll a 


-0.10 


-0.16 


-0.47 


Sample day Ul 


0.02 


-0.02 


-0.18 


Monthly Ul 


-0.26 


-0.61 


-0.23 


Davidson season 


0.06 


0.19 


0.15 


Oceanic season 


0.31 


0.41 


0.03 


Upwelling season 


-0.21 


-0.51 


-0.10 


Year of study 


-0.14 


-0.35 


-0.09 


MBAY Vib7 


0.01 


0.01 




MBAY Vib4 


0.55 










Total Vibrio 


MBAY Vib4 


MBAY Vib7 


All phytoplankton 


0.18 


-0.44 


-0.38 


Total diatoms 


0.10 


-0.53 


-0.43 


Total dinoflagellates 


0.17 


-0.17 


-0.43 


Synechococcus 


0.02 


0.20 


0.41 


Surface temp. 


0.21 


0.45 


0.07 


Surface Chlorophyll a 


-0.21 


-0.30 


-0.29 


Surface nitrate 


-0.20 


-0.39 


0.00 



Significant parameters (p < 0.05) are in bold type. Ul, upwelling index. 
(A) Upwellling parameters verified against surface data from the M2 mooring 
station. (B) CTD data and environmental parameters. (C) Surface phytoplankton 
and M2 mooring station data. 



within 35 defined regions, 28 of which are located between the 
primers used to construct the clone library for this study (Wiik 
et al., 1995; Jensen et al., 2009). It is thus not surprising to 
have identified a significant number of OTUs and to observe 
deep branching within the phylogenetic tree of the MB clone 
library. Intraspecies sequence variability can also be high within 
Vibrio species (Jensen et al., 2009) and the identification of Vibrio 
splendidus within different branches of the tree is supported by 
observations of up to 2% difference within the 16S rRNA gene of 
this species (Jensen et al, 2009; Le Roux et al, 2009). Defining 
relationships within Vibrio populations is even more complex 
as their genomes contain multiple copies of the rRNA operon 
(Heidelberg et al, 2000). Thirteen copies have been identified 
in V. natriegens and 12 copies in both Allivibrio fischeri and A. 
salmonicida (Lee et al, 2009). Intragenomic heterogeneity within 
these copies may be high, as seen in V. parahemolyticus (Harth 
et al., 2007), or non-existent as is noted in other Vibrio species 



(Coenye and Vandamme, 2003). Since this study was not designed 
to differentiate between heterogeneous intragenomic copies and 
individuals, our results may overestimate diversity within the MB 
clone library. With that being said, sequences did cluster with over 
20 known Vibrio species suggesting that MB does in fact contain 
diverse Vibrio populations. 

Multiple copies of the rRNA operon may also affect the 
qPCR estimates of abundance presented in this study, but is 
not expected to significantly change the interpretations of our 
findings. Even an estimate of 20 gene copies per cell results in esti- 
mates of abundances that fall within the range of concentrations 
observed in other coastal waters (Thompson et al, 2004a), and 
are higher than observations from samples in similar tempera- 
ture ranges (Randa et al., 2004; Eiler et al., 2006). Temperature 
is often cited as a significant driver in Vibrio population abun- 
dance (Heidelberg et al., 2002; Maeda et al, 2003; Eiler et al, 
2006, 2007; Hsieh et al, 2008), and while temperature variabil- 
ity was correlated to the two subpopulations examined in this 
study, it was the upwelling index (MBAY Vib4) and depth and sta- 
tion (MBAY Vib7) that were more significant factors in defining 
subpopulation abundance. 

MBAY Vib4 and MBAY Vib7 are separate subpopulations 
of the Vibrio community with distinct niches in the MB. The 
MBAY Vib4 is defined by close association with Allivibrio sp. 
and represents an oceanic season subpopulation. This season is 
characterized by a lack of upwelling, lower sea surface salini- 
ties, warmer sea surface temperatures and waters that are influ- 
enced by wind relaxation events and the slow flowing California 
Current (Rosenfeld et al., 1994). MBAY Vib7 represents a deeper 
water subpopulation of Vibrio sp. with a greater association with 
the offshore waters at station M2. This station is influenced 
by upwelled waters mixed with California Current waters, the 
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FIGURE 3 | Neighbor-joining phylogenetic trees of Vibro sequences 
from Monterey Bay at stations C1 and M2. Tree constructed in ARB from 
partial 16S rRNA sequences (£. coli 27-680 bp) using ARB neighbor-joining 

(Continued! 



FIGURE 3 | Continued 

distance matrix method with a Jukes-Cantor correction. T, total number of 
sequences in clade, MB, MB sequences in clade. Sample number, station, 
depth, and date are listed for Monterey Bay sequences. (A) Major groups 
of environmental Vibrio sequences recovered from the MB. Shewanella 
algae included as the tree root. (B) Allivibirio and MBAY Vib4 group 
sequences, and (C) Vibrio penaecida and MBAY Vib7 sequences. Scales 
represent 10% or 0.1 % nucleotide substitutions. 



California Undercurrent and the Davidson Current-a northward 
flowing current that develops in winter months (Tisch et al., 
1992; Breaker and Broenkow, 1994; Rosenfeld et al, 1994). The 
California Undercurrent usually flows below 100 m along the 
California shelf but often nears the surface during the Davidson 
season (Tisch et al, 1992; Pierce et al, 2000; Tseng et al, 2005). 
The largest hydrographic changes around MB occur during tran- 
sitions between upwelling regimes (i.e., during oceanic phases) 
and in winter when the horizontal and vertical thermal gradients 
are reduced due to the northward flow of the Davidson Current 
(Bac et al, 2003; Storlazzi et al, 2003; Warn-Varnas et al., 2007). 
Dynamics of the MB Vibrio populations examined in this study 
seem to reflect these hydrographic processes, and may be highly 
influenced by them-MBAY Vib4 during the oceanic phases and 
MBAY Vib7 during the winter phases. 

In 2006, the MB had experienced 4 years of anomalous oceano- 
graphic conditions. Delayed and unusually shallow upwelling 
affected the food web in MB up through higher trophic levels 
(Goericke et al, 2007). Changes included a shift in dominant 
toxin-producing algal species from diatoms to dinoflagellates, 
poor recruitment of krill and low zooplankton biomass as well 
as seabird reproductive failure (Peterson et al, 2006; Goericke 
et al., 2007; Jester et al, 2009). MBAY Vib4 was positively cor- 
related to this year possibly due to delayed upwelling, reduced 
chlorophyll a concentrations and warmer surface temperatures 
(Table 3). Increased abundances of Vibrio populations could also 
have been influenced by decreased competition from phytoplank- 
ton for resources or reduced top down control from zooplankton. 
Following the 2002-2006 lull events, unusually strong upwelling 
was observed in 2007 (Kaplan et al., 2009), and in 2008 a strong 
development of upwelling was observed at the beginning of the 
season (Figure 3). MBAY Vib4 populations showed significantly 
reduced abundance during those upwelling seasons. MBAY Vib7 
and the total Vibrio populations did not seem to be strongly 
affected by these events. 

Overall, the Vibrio populations examined in MB seemed to be 
influenced by larger scale upwelling events and shifts in currents 
and oceanographic seasons rather than individual environmental 
factors. In the future, more extensive analysis of nutrient availabil- 
ity and oceanographic parameters (i.e., flow velocities and wind 
patterns) may provide better insight into how upwelling charac- 
teristics and water mass flows play a role in Vibrio population 
dynamics. Since Vibrio populations might have significant influ- 
ence over nutrient availability through recycling chitin and other 
bound nutrient resources, understanding the role of these bac- 
terial populations in both surface and deeper waters can better 
our understanding of the productivity potential within upwelling 
regions. 
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FIGURE 4 | Monthly Upwelling Index (m 3 s- 1 100 m of coastline- 1 ), temperature (°C), salinity (ppt), chlorophyll a (u.g m 3 ) and total Vibrio, MBAY Vib4 
and MBAY Vib7 concentrations (natural log of gene copies liter -1 ) for stations CI (A) and M2 (B). Black dots represent discrete samples within the CTD cast. 
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